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1 . SUMMARY 


The photoelectrochemical effects that occur at an illuminated semiconductor- 
liquid interface may provide a unique means of enerqy conversion and storage. A 
semiconducting electrode illuminated by light with energy greater than the semicon- 
ductor bandgap energy may induce a chemical reaction among the constituents of the 
interfacial region. A model of the charge transport processes in the semiconductor 
is presented in this report. This model is based upon semiconductor device theory, 
and it incorporates, in an essential way, the nonlinear processes characterizing the 
diffusion and reaction of charge carriers in the semiconductor. The model is used to 
study conditions limiting useful energy conversion, specifically the saturation of 
current flow due to high light intensity. Numerical results describing charge dis- 
tributions in the semiconductor and its effects on the electrolyte are obtained. 

Results of experimental efforts are also summarized. These results include the 
following: 

1 . An estimate based on experiments of the maximum rate at which a semiconductor 

photoelectrode is capable of converting electromagnetic energy into 
chemical energy 

2. The effect of cell temperature on the efficiency of a photoelectrochemical 

cell 

3. A method for determining the point of zero zeta potential (pzzp) for 

macroscopic semiconductor samples 

4. A new technique using platinized titanium dioxide powders and ultraviolet 

radiation to produce chlorine, bromine, and iodine from solutions 
containing their respective ions 

5. The results of experiments characterizing a new class of layered compounds 

called transition-metal thiophosphates for use as photoelectrodes in 
photoelectrochemical cells 

6. A discussion of a technique used to produce high conversion efficiency from 

laser radiation to chemical energy. 

2. INTRODUCTION 

In 1972, Fujishima and Honda (refs. 1 and 2) demonstrated that sunlight, inci- 
dent upon an n-TiC^ electrode, would catalyze the electrolysis of water. This 
research stimulated much more study (refs. 3 through 7) of the photoelectrochemical 
processes occuring at a semiconductor-electrolyte interface and led to the discovery 
of a variety of useful chemical reactions stimulated at the surface of a photosensi- 
tive electrode. There is a growing body of evidence to suggest the possibility of 
developing a semiconductor electrode device capable of converting radiant energy into 
electrical energy or into a chemical fuel. If such a device could operate effi- 
ciently and on a large scale, there would be important applications both on Earth and 
in space. Specifically, there may be some unique applications of these devices in 





the environment of space, where they can be used as a power supply or as part of a 
propulsion system. One particularly interesting concept is that hydrogen and oxygen 
produced from solar radiation by a semiconducting electrode cell may be used to 
replenish the stores of these chemicals consumed in a fuel cell. This would provide 
a closed-cycle solar-energy system. It is also interesting to conceive of this type 
of cell as a receptor of energy transmitted from one point in space to another by a 
laser beam as part of a space energy -distribution network. 

Figure 1 is a representation of a semiconductor-electrode electrochemical cell. 
The critical element in such a cell is the combination of the photosensitive semicon- 
ductor electrode and the semiconductor-electrolyte interfacial region. It is within 
this region that the incident photon energy is transferred throuqh various energy 
levels in the semiconductor and in the electrolyte, until it eventually appears in a 
useable form as chemical or electrical energy. The overall process of converting the 
incident radiant energy to its final form is the result of a series of simpler pro- 
cesses of a physical or chemical nature. In order to understand the operation of a 
complete photoelectrochemical cell, it is necessary to understand the individual 
processes which occur, as well as the way each of the separate processes affect each 
other. In this report, a theoretical treatment of the relevant processes taking 
place within the semiconductor electrode and at its interface with the electrolyte is 
described. This treatment is based upon a rigorous application of the macroscopic 
theory of charge generation, recombination, and transport in semiconductors. The 
equations which are derived for charge transport in a semiconductor electrode are 
nonlinear. Therefore, an exact solution was not obtained. A systematic procedure 
for obtaining approximate solutions to them are developed, and numerical results 
based upon the approximations are presented. The previous theoretical work was 
developed from a linearized theory. There is an inherent limitation in this 
approach, in that results from the linearized theory do not apply in all situations. 
There are some circumstances where the nonlinearity of the processes is important, 
and one can hope to glean some information about these circumstances from the more 
rigorous development presented herein. Treatment of the saturation of photoeffects 
in the semiconductor at high light intensity is one illustration of this point. An 
alternative method for discussing the transport of excess charge carriers is pre- 
sented in references 7 and 8. It should also be noted that although the specific 
case of a semiconductor electrode is addressed, the theory developed herein is appli- 
cable, with minor modifications, to other semiconductor devices. 

The major regions of interest within an electrochemical cell are shown in the 
schematic in figure 1. The phenomena responsible for photos timulated catalysis or 
for the production of electricity in a semiconductor-electrode cell are controlled by 
processes taking place in the semiconductor-electrolyte interfacial region. A 
detailed and thorough understanding of all these interrelated processes is lacking at 
present, but the following description provides a sketch of the probable events 
occurring in the interfacial region. When a semiconductor electrode is brought into 
contact with an electrolyte, there is a redistribution of the mobile charge within 
the semiconductor and of the ionic charge within the electrolyte. These redistribu- 
tions result in a steady-state space-charge distribution throughout the interfacial 
region. The space-charge distribution within the electrolyte interfacial region has 
been studied extensively (refs. 9 through 12) and can be described adequately by the 
Gouy-Chapman-Stern theory of ionic-charge distribution at solid-liquid interfaces. 

The extent of this charge distribution and the net charge contained in it determine 
the net charge on the semiconductor. The space-charge distribution within the semi- 
conductor is the source of an electric field near the surface, and this field plays 
an important role in the photosensitivity of the cell. Illumination of the semicon- 
ductor by light of sufficiently high energy produces electron-hole pairs in the semi- 
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conductor. These pairs are then separated by this space-charge field. Part of the 
photoproduced charge is drawn to the surface of the semiconductor, where it can 
interact (via electron transfer) with the electrolyte. This interaction causes a 
chemical reaction among the constituents of the interfacial region. 

The effect of the redistribution of charge in the semiconductor electrode when 
it is placed in contact with the solution is to produce an electric field near the 
surface. This field gradient alters the semiconductor energy bands near the surface 
(ref. 13). If the electrode and the electrolyte are placed in contact, the energy 
levels are as shown in figure 2(a). This figure illustrates the case of an n-type 
semiconductor electrode. The Fermi energy of the electrolyte Ep RE DOX (chemical 
potential) is determined by the solvent and the ionic species in it. When the semi- 
conductor electrode and the electrolyte are placed in contact, there is a readjust- 
ment in charge which brings the two Fermi energy levels together at a value extremely 
close to Ep redox* surface values of the conduction band edge E c and the 

valence band' edge E v remain fixed, and the final energy -band picture is shown in 
figure 2(b). An important special value of e f, REDOX that for w hich the band 

bending is zero, the flat band energy Ep r EDOX * ' J ^ ie ^tference between Ep REDO x 
and Ep RE DOX amount °f band bending induced by the electrolyte. By using 

this band-bending picture, the operation of a photoelectrochemical cell can be repre- 
sented as in figure 3. A photon, absorbed by the semiconductor, promotes an electron 
from the valence band to the conduction band and leaves a hole in the valence band. 
The photoproduced charges are separated by the semiconductor space-charge field. The 
electron moves through the bulk of the semiconductor to the external circuit, and the 
hole moves to the surface, where it can interact with the electrolyte in an oxidation 
reaction. The current due to the electron motion through the circuit is compensated 
by an ionic current in the electrolyte to the counter electrode, where it partici- 
pates in a reduction reaction. 

The electronic states of ions in solution can be represented as electron energy 
levels, subject to fluctuations in time due to the thermal motions of polar molecules 
in the solution (ref. 14). Thus, there are at least two broadened electronic levels 
in the electrolyte corresponding to redox reactions within the cell. There may be 
additional energy levels corresponding to electronic states within the semiconductor. 
One level represents the energy of dissociation, E^, and the other levels may repre- 
sent possible surface states. All these energy levels in the interfacial region are 
shown in figure 4. Although some semiconductor-electrode electrochemical cells may 
produce relatively large photocurrents, others decompose irreversibly because of the 
redox reactions at their surface. This decomposition occurs when the energy of dis- 
sociation E_ is more negative than E_ so that the dissolution of the elec- 

trode is energetically possible. The electrode shown xn figure 4 will not dissolve. 
There are certain semiconducting metal-oxide photoelectrodes that are inherently 
stable (refg. 3 through 5). Unfortunately, some of the stable semiconducting photo- 
anode materials (Ti0 2 and SrT i° 3 / e.g.) have bandgaps so large that only ultraviolet 
radiation ( »3 percent of the solar spectrum) is absorbed. However, the energy con- 
version efficiency (photo-to-chemical) is quite high. It has been found recently 
that an appropriate electrolyte often stabilizes semiconducting electrodes that were 
thought to be unstable (ref. 3). Consequently, electrodes with better response to 
the solar spectrum might be developed. These stabilized semiconductor-electrolyte 
junctions have not yet stimulated electrolysis. 
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3. SYMBOLS 


n 

°P 

E 

EA 


D A 

* 


E, 


Ed 

* 


T,I 


J F, REDOX 


Ph 


^SS 

Ev 


n _ 2 

cross-sectional area of electrode, m^; also a matrix of constants, s 

dimensionless parameter characterizing bulk concentrations of electrons and 
holes in semiconductors; also an arbitrary parameter 

numerical index 

speed of light, m-s _1 

matrix with diffusion coefficients as elements, m -s 

9 — 1 

diffusion coefficient for electrons, m -s 

2 — 1 

diffusion coefficient for holes, m -s 
energy, eV 

electron affinity, eV 

energy level of an acceptor atom, eV 

effective ionization energy of an acceptor atom referred to vacuum 
energy, eV 

lowest energy in conduction band, eV 
energy of dissociation, eV 

effective ionization energy of a donor atom referred to vacuum energy, eV 

Fermi energy in a semiconductor, eV 

Fermi energy of an intrinsic semiconductor, eV 

Fermi energy in electrolyte at flatband condition, eV 
bandgap energy, eV 
energy of bent bands, eV 
effective energy of bent bands, eV 
photon energy, eV 

energy level of single electronic state, eV; also magnitude of electric 
field strength at semiconductor surface, V-nT 1 

effective energy level of single electronic state, eV 

energy of surface states, eV 

highest energy in conduction band, eV 
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( ) 


E(k) 
E(x) 
e 

FF 
P n 

F P 
F(u) 

P 1/2 
F( ) 

G n 

G P 

g 

g* 

g-r 

1 o 

H( , ) 

h 

I 

K ) 

X n 

X P 

I ph 

X 0 

X 0,s 


U 


electronic energy, eV 
electric field strength, V-m -1 
electronic charge, C 
fill factor, dimensionless 

_ O _ 1 

net source of electrons, m -s 
net source of holes, nT^-s" 1 

matrix with elements giving net source of charge carriers in a 
semiconductor, itT^-s -1 

Fermi integral, dimensionless 

Fermi distribution, dimensionless 

_ O _ 1 

net rate of generation of electrons, m -s 

— 3 1 

net rate of generation of holes, m -s 

number of degenerate states of a vacant electronic level 
number of degenerate states of an occupied electronic level 
rate of photoinduced generation of charge carriers, m”^-s”"' 
rate of thermal generation of charge carriers, m""^-s~1 
first integral of Poisson-Boltzmann equation, dimensionless 
Planck constant, eV-s 
electric current, A 

matrix giving inhomogeneous part of transport equation, C-s-m” 2 
inhomogeneous term in transport equation -tor electrons, C-m”^ 
inhomogeneous term in transport equation for holes, C-m”^ 
photocurrent, A 

intensity of incident light, m“ 2 -s -1 

— 9 i 

saturation level of incident light intensity, m -s 

— 2 ^-1 

electron current density, C-m -s 

— 2 —1 

hole current density, C-m -s 

current density due to absorption in the Urbach region, C-m“ 2 -s -1 
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n > 

* 

m 


m 


n 
* 

m P 

NHE 

N A 


N_ 


N, 


n i 

P 

Pa 

Pb 

Q ± 

q 

R 

R n 

R P 

SC 

SCE 


Boltzmann constant, J-K ^ ; also particle wave number, m ^ 

Debye length, m 

electron diffusion length, m 

hole diffusion length, m 

effective Debye length, m 

linear functional 

effective mass, kg 

electron effective mass, kg 

effective hole mass, kg 

normal hydrogen electrode 

concentration of acceptor atoms, m -3 

conduction band concentration factor, m -3 

concentration of donor atoms, m *"" 3 

valence band concentration factor, m ~ 3 

electron concentration, m ~ 3 

bulk value of electron concentration, m “ 3 

concentration of electrons in donor atoms, m 

, . . — ^ 
intrinsic concentration, m 

hole concentration, m “ 3 

concentration of holes in acceptor atoms, m “ 3 
bulk value of hole concentration, m ~ 3 

_ n 

surface charge excess, C-m 
electronic charge, C 

— 1 — 1 

gas constant, J-kmol -K 

rate of recombination of electrons, m"’ 3 -s ” 1 
rate of recombination of holes, m~ 3 -s _1 
semiconductor 

saturated calomel electrode 
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T 

t 

u 


V U 1 


v 


app 


V, 


bb 


V 


ext 


V, 


FB 


V 


V 


op 


V 


prod 


V 


V 
W 

x,x 


0 


y 

yi/2 


temperature , K 
time, s 

normalized electron energy, dimensionless; also matrix with charge carrier 
concentration as elements, m “ 3 

bulk value of normalized electron energy, dimensionless 

• — 

• iterated solutions to transport equation, m 

potential difference, V 

applied potential difference, V 

potential of bent bands, V 

potential supplied by external source, V 

potential of the Fermi level, V 

flatband potential, V 

potential across Helmholtz layer, V 

electron velocity, m-s -1 

onset potential, V 

hole velocity, m-s“^ 

difference in redox potential of half -cell reactions, V 
surface potential, eV 

surface potential, V; also redox level for hydrogen, V 

bulk potential, eV 

depletion width, m 

dimensionless parameters 

distance within semiconductor, m 

penetration parameter, dimensionless 

special value of penetration parameter, dimensionless 
absorption length for incident radiation, 
electron drift parameter, m^ 
hole drift parameter, m 
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“o 

p 

A 

AE 

A „ 

AG 

AT 

AV 

4v pb 


c,,e 2 


^nax 

X 


P C (E) 

P V (E) 


Urback absorption length, m ” 1 
thermal energy factor, C-J - ^ 

normalized electron potential, dimensionless 
change in energy, J 

position of Fermi level with respect to the conduction band edge, J 

change in Gibbs free energy, J 

change in temperature, K 

change in potential, V 

change in flatband potential, V 

normalized electron energy, dimensionless 

electrostatic permativity, v - 1 -m -1 

numerical parameters 

electrochemical cell efficiency, dimensionless 

maximum value of electrochemical cell efficiency, dimensionless 

wavelength, m; also a dimensionless parameter 

chemical potential, eV; also drift mobility, m 2 -V~ 1 -s -1 

electron drift mobility, m -V -s 

hole drift mobility, 

photon frequency, Hz 

charge distribution, C - m 

density of electronic states in the conduction band, m” 3 -eV ” 1 

density of electronic states in the valence band, m” 3 -ev “ 1 

electron lifetime, s 

hole lifetime, s 

electrostatic potential, V 

bulk value of electrostatic potential, V 
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X 


electron affinity, eV 


u> n electron transport parameter, dimensionless 

0)p hole transport parameter, dimensionless 

Prime indicates derivative with respect to x. 

An arrow over a symbol indicates a vector quantity. 


4. SPACE-CHARGE DISTRIBUTION IN A SEMICONDUCTOR ELECTRODE 
4.1 General Description of Electrode 

In this section, a foundation is laid for a semiclassical analysis of the phys- 
ical processes occurring in an illuminated semiconductor electrode. The basic system 
under investigation is depicted in figure 5, which schematically represents the 
semiconductor -electrolyte interfacial region. The presence of ionic charge in the 
electrolyte and possibly a distribution of charge on the semiconductor surface cause 
the mobile charge within the semiconductor interfacial region to be displaced, leav- 
ing a distribution of fixed ionic charge. These distributions of fixed and mobile 
charges near the semiconductor surface are the source of an electrostatic field. In 
the sections which follow, we examine the form of this field and subsequently study 
the role of this field in the photoeffects occurring at the junction. The assump- 
tions which underlie our analysis are as follows: 

1 . The semiconductor is a single-crystal cylinder with one base plane of surface 
area A in contact with the electrolyte and with the opposite base plane in ohmic 
contact with the metallic counter electrode of the cell. The diameter of the crystal 
base is much larger than the crystal thickness. Thus, the variables that charac- 
terize the electrical properties of the crystal depend upon a single spatial param- 
eter, the depth from the semiconductor-liquid junction. An n-type semiconductor is 
assumed; this is most frequently the case in practice. It is a simple matter to 
adjust the results to apply to a p-type material. 

2. The doping of the semiconductor is attributed to a concentration of acceptor 

and donor ions and distributed uniformly throughout the crystal. Further- 

more, the separation between acceptor and donor ions is large enough to be able to 
ignore interaction between ions. On the other hand, the spacing between acceptor and 
donor ions is small enough, compared with the space -charge penetration depth, that 
the discrete ionic-charge distribution can be represented as a continuous 
distribution. 

3. If a surface state is charged, it resides on the surface and not within the 
semiconductor. This implies that the space-charge distribution is caused by fixed 
ionic charge and mobile charge carriers only. 

4. The semiconductor is nondegenerate; that is, the Fermi energy differs from 
the conduction and valence band energies by at least 3kT. This assumption is dis- 
cussed more fully in the appendix. The most prominent effect of this assumption is 
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that the concentration of electrons and holes can be described by a Boltzmann distri- 
bution. Namely, the concentration of electrons is given by 


n 



(4.1 ) 


where 


N 


C “ 



3/2 


and the concentration of holes is given by 


P 


N y exp 



kT 


(4.2) 


where 



If the semiconductor is maintained in thermal equilibrium, the Fermi energy 
E f maintains a constant value at all points. Consequently, if the concentration of 
carriers n and p changes from point to point in the semiconductor as it does near 
the semiconductor electrolyte junction, this variation must mean a local change in 
E c and Ey . That is, the energy bands are bent. This bending is represented sche*^ 
matically in figure 6, which depicts the band bending of an n-type semiconducting 
electrode. Ihe band bending is a result of various mobile charges in the junction 
region coming into electrostatic equilibrium. This charge distribution is described 
in the following section. 


4.2 Derivation of Conditions of Electrostatic Equilibrium - 
Poiss on-Bo ltzmann Equation 

Let the uniform concentration of acceptor and donor atoms be represented by 
N a and N d , respectively. Assume that these atoms are completely ionized, a circum- 
stance readily satisfied at room temperature. (See appendix.) Then, if n(x) and 
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p(x) represent the carrier concentrations as a function of depth in the semiconduc- 
tor, the net charge distribution within the semiconductor is given by 


p(x) = e [p ( x ) - n (x) + N d - N^] 


(4.3) 


Deep within the semiconductor, the charge carrier distributions achieve their 
bulk values, and charge neutrality obtains. This limiting behavior is expressed by 


lim 

x-ko 


n(x) 


n 


b 


lim 

x-**> 


p(x) 


= Pv 


and 


lim 


p(x) 


0 


Combining the results, we find that the charge neutrality condition away from 
the junction requires that 


Pb - n b + n d - N A = 0 


(4.4) 


Using this last result, we can express the net charge distribution in terms of 
the carrier excess over its bulk value by 


p(x) = e{[p(x) - P b ] - [n (x) - n b ]} 


(4.5) 


This charge distribution is taken to be the source of the equilibrium electrostatic 
field within the semiconductor. 

Whenever the Fermi energy level within the semiconductor differs by more than 
3kT from either band edge, the electron and hole concentrations at any point must 
conform to the Boltzmann distribution (eqs. (4.1) and (4.2)). By multiplying the two 
distributions together, we obtain the equation 


np 


N c N v 



(4.6) 


Since the right-hand side of this expression is constant (for a fixed temperature), 
the product of n(x) and p(x) must also be constant. Equation (4.6) can be used 
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to introduce a convenient normalization for these concentrations. In the special 
case of the bulk region of an intrinsic semiconductor where there is no ionic charge , 
the condition of charge neutrality requires that electron and hole concentrations be 
equal, and the concentration in this case n^ is used as a norm. Since equa- 
tion (4.6) holds in all cases (including this special one), n^ is given by 


n 


i 


2 


N c N v 



(4.7) 


Using n^ as defined above, the concentration of electrons and holes in the general 
case can be normalized as 


n 


n. 

l 


exp 



- E*(x) 
kT ~~ 


(4.8) 


and 


p = n^ exp 


E p - E I (x) 

icT 


(4.9) 


where 


* 



= c + V 



(4.10) 


The relation between and the bent energy bands is shown in figure 6. The elec- 

trostatic energy of an electron within the semiconductor relative to an electron in 
the intrinsic bulk is given by 


e <J)(x) = E p - E l (x) 


(4.11) 
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In terms of this electrostatic energy , the expressions for the carrier concentration 
become 


n(x) = n. exp[^j^-] 


and 


p(x) 


= n ± exp jj- 


e <t> ( x)~ | 

kT J 


( 4 . 12 ) 


( 4 . 13 ) 


Now the limiting value of the electrostatic potential within the semiconductor 

defines the bulk value $ as 

b 


lim 

x-x» 


<t>(x) 



E*(co) 


Comparing this value to the bulk values of the concentrations yields 


n 


b 


n . 
r 



( 4 . 14 ) 


and 


P 


b 


n . 
l 



( 4 . 15 ) 


The electron and hole concentrations can also be expressed in terms of the bulk 
values by means of the equations 



( 4 . 16 ) 


( 4 . 17 ) 


1 3 



By using the expressions for n(x), p(x), n bf and p fe , normalized by the 

intrinsic concentration n^, the concise expression is obtained for the charge dis- 
tribution in terms of the electrostatic potential, namely 


p( <j>) = 2en^ 




( 4 . 18 ) 


The condition for electrostatic equilibrium (Poisson equation) becomes in this case 
the second-order, nonlinear differential equation 



( 4 . 19 ) 


sometimes referred to as the Poisson-Boltzmann equation. By means of the following 
changes in variables 


u 


kT 


and 


u 


b 


e 


*b 


kT 


this equation is written in the form 


u 


H 


— —(sinh u - sinh u b ) 
L 

D 


( 4 . 20 ) 


/here is the Debye length given by the relation 


0 2 
2e n . 

i_ 

skT 


( 4 . 21 ) 


14 



4.3 Solutions to Poi s son-Bo ltzmann Equation 

This section contains a description of the properties of solutions to the non- 
linear differential equation 


u"(x) 


— {sinh[u(x)] 
L D 


sinh ( u b )} 


(4.22) 


defined for 0 < x < » and satisfying the limiting conditions 


u ( 0) = u Q 


and 


lim 

x-*» 


[u(x) 


“J ■ 0 


To examine the asymptotic behavior of solutions to equation (4.22), the function 6, 
as defined by the equation 6 = u(x) - u b , is introduced. The equation satisfies the 
differential equation 


cosh(u ) sinh(u ) 

6" = — sinh( 6) + [cosh(6) - l] 


For small values of 6, this equation assumes the linear form 


6 " = 


cosh(u ) 


and has a bounded exponential solution with a characteristic asymptotic length given 
by 



(4.23) 
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Thus, bounded solutions to equation (4.22) display an asymptotic behavior of the form 


u(x) -* u^ + A exp^- 

A first integral (energy integral) to the Po is s on-Bo ltzmann equation can be obtained 
directly and leads to the first-order equation of the form 


[6’ (x) ] 2 = — {cosh(6) 

L 2 

CD 


+ a[sinh(6) - 6] } 


(4.24) 


where 


a = tanh(u ) = 


n b - P b 
n b + P b 


(4.25) 


The right-hand side of equation (4.24) is positive for x > 0, and it has a single 
zero at x = 0; the root of that equation thus yields the equation 


6' (x) = ± H[6(x),a] 

GO 


The function H(5,a) is given by the expression 

H(6,a) = (2{cosh(6) - 1 + a[sinh(6) - S]}) 1//2 


(4.26) 


(4.27) 


The ambiguity in sign resulting from the square-root function can be resolved as 
follows. Since H(6,a) is positive or zero for values of 6 and a within their 
appropriate range of variation, the sign in equation (4.26) determines the slope of 
6, and this slope maintains the same sign for all values of 6. Furthermore, the 
value of 6 must approach zero within the bulk. Thus, in the case of the plus sign, 
6 approaches zero monotonically with increasing y from a negative surface value, 
as occurs in a p-type material depletion region. In the case of the minus sign, 6 
approaches zero monotonically with increasing y from a positive surface value, as 
occurs in an n-type material depletion region. 

By making a scale change to the dimensionless independent variable y = x/L , so 
that 6(x) becomes 6(L y) = A(y), one of the two remaining parameters in this equa- 
tion can be eliminated, and the following equation can be obtained. 


A' (y) = ±H[A(y) ,a] 


(4.28) 
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This equation depends upon the single parameter a. Integrating equation (4.28) 
yields 



ds 

H(s f a) 


±y 


(4.29) 


This integral gives the value of A as an implicit function of the penetration 
parameter y in terms of the two parameters A(0), which specifies the potential 
drop across the semiconductor normalized to mean thermal energy 


A(0) =f^[(t>(0) - *(«)] 


(4.30) 


and a, which specifies the doping fraction and is defined in equation (4.25). The 
penetration parameter gives the actual penetration depth in units of the asymptotic 
length given in equations (4.21) and (4.23). 

Results of a numerical evaluation of equation (4.29) are shown in figure 7. The 
dimensionless penetration parameter is plotted along the abscissa. No origin of this 

coordinate is defined in the graph , so that the zero value of y can be established 

beneath any specified initial value of A(0). The values of A(y) are plotted along 

the ordinate in a logarithmic scale increasing away from a central zero value in both 

directions. Naturally, the actual value of zero does not appear in the graph. 

Curves shown on the upper half of the figure are associated with the minus sign in 
equation (4.28) and represent a depletion region? curves in the lower half are asso- 
ciated with the plus sign and represent an accumulation region. 

Under special circumstances, the solution to the Poisson-Boltzmann equation 
given by equation (4.29) takes on a simple form. For example, if the parameter a 
is equal to zero, the integral in equation (4.29) can be evaluated analytically to 
obtain the expression 


tanhj^— = exp(±y) tanhj ^ — (4.31) 

When A is either extremely small or extremely large, approximate values of the 
integrand can be developed which can be evaluated analytically. If, for example, 

|x| « 1, H(x,a) can be approximated by 


H(x,a) « x 



and the resulting integral in equation (4*29) can be evaluated over intervals where 
|a| < 1 to yield 



then for small initial values of A(0), A(y) changes exponentially to yield 


A(y) * exp(±y) 


Thus, all the curves in figure 7 approach a straight line with a slope of ±1 
as y decreases to zero. Similarly, if |x| >5, the value of H(x,a) becomes 


H(x,a) 


» 


( 1 + a ) 


1/2 



Integrating this function over an interval where 


| A | >5 yields 



±(1 


+ a) 


1/2 


y 

2 


This explains the steep descent or ascent of the curves in figure 7 for large values 
of A. These curves are descending or ascending from an initial value A(0) with a 
half-length of 



4.4 Abrupt Approximation 


The solutions to the Poisson-Boltzmann equation indicate that the electrostatic 
potential in the space-charge region of the semiconductor decreases monotonically (or 
increases monotonically in the case of an accumulation layer) from its surface value 
to its value in the bulk over distances comparable to L . For some analyses, it is 

OO 

convenient to approximate the electrostatic properties of the space-charge region 
still further by assuming that the electrostatic field changes linearly from its 
value at the surface to zero at a fixed depth W within the semiconductor and that 
the field remains zero throughout the bulk. In this abrupt approximation, the field 
strength is nonzero only within the region of width W, called the depletion width, 
at the surface. (See fig. 8.) Again, the depletion region of an n-type semiconduc- 
tor is considered. The source of the space-charge field is the ionic charge distri- 
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bution p, exposed by the depleted majority carriers. From the assumed field distri- 
bution (fig. 8) and from the Poisson equation, it can be seen that 


E 

dE s_ _ _ 

dx W £ 


From this, the shape of the charge distribution is obtained as 


P = 



(4.32) 


for 0 < x < W and p = 0 elsewhere. The potential drop across the depletion layer 
is then given by 


AV = 



(4.33) 


Eliminating E g from the two previous equations shows that the depletion width is 
determined by the potential drop across the depletion layer 


2 2eAV 

W = 

P 


(4.34) 


If it is further assumed that all donor atoms are ionized (a reasonable assumption 
for temperatures near 300 K), then 


p = eN = en 
H D b 


Finally, the following expression relates the depletion width W to the potential 
drop AV: 


W = W q (AV) 1//2 (4.35) 

where 
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It follows that the mean electric field in the depletion region is given by 


E 



(AV) 


1/2 


W„ 


(4.36) 


In this abrupt approximation, the mean space-charge field varies as the square root 
of the potential drop across the region, and it penetrates to a depth Wq when the 
potential drop equals 1 V. 

In some of the analyses which follow, it is expedient to use the abrupt approxi- 
mation for the space-charge field of a depletion region instead of a field derived 
from the more general conditions of electrostatic equilibrium. It is of interest to 
compare the measure of field width in the abrupt approximation W Q with the cor- 
responding measure of the Poisson- Boltzmann field L^. As shown previously, 
is given by 


L 


2 


CO 


£ ^ T (cosh 

2e n. 

3 . 



-1 


where, from equation (4.14) 


n b + p b 

COSh U b = 2n". 


In the depletion region of an n-type semiconductor n^ » p fo ; therefore, is rela- 

tive to Wq through the relation 



At room temperature, this relation implies that W Q » 9L OT * 


(4.37) 


5. CHARGE TRANSPORT IN A SEMICONDUCTOR ELECTRODE 
5.1 Transport Equations 

As shown in the preceding section, there are three interdependent variables 
involved in the macroscopic description of current flow in a semiconductor. Two of 
these are the concentrations of mobile charge carriers n(x,t) and p(x,t), repre- 
senting the local concentration of electrons and holes, respectively. The third is 
the local electric-field intensity E(x,t). The electric field within the semicon- 
ductor is composed of two parts - the external field imposed from outside the system, 
and the internal field whose source is the internal charges (both fixed and mobile) 
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within the semiconductor-liquid interfacial region. These three quantities vary most 
in the interfacial region, and within the bulk of the semiconductor assume constant 
limiting values. The bulk electric-field intensity is taken to be zero, as its value 
is small compared with field values near the surface. The bulk values for the elec- 
tron and hole concentrations are given by n^ and p fa , respectively. 

The free charges move under the influence of local electric fields (drift) and 

local concentration gradients (diffusion), so that the current density for each car- 
rier is given by 

J = efnV + D Vn) (5.1 ) 

n v n n ' 

and 

J = efpV - D Vp) (5.2) 

p ^ p P 

In these expressions, D n and are diffusion coefficients for the indicated 

carriers, and V n and are corresponding drift velocities. Throughout this 

discussion, it is assumed^ that the field intensity is small enough to allow the 
linear approximation given by 


-► 

V 

n 



and 


-y 

V = 
P 


-► 
LL E 

JP 


(5.3) 


(5.4) 


The carrier mobilities \x and (j, are thus related to the diffusion coefficients 
through the Einstein relations ^ 


D 

n 



(5.5) 


where 


^ kT 

The continuity relations for the carrier transport are 


5n J_ 
at e 


-y 

J 

n 


G - R 


n n 


(5.6) 
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and 


5n 

at 



} 

P 


G 

P 


R 

P 


(5.7) 


The terms on the right-hand side of each of these formulae represent the net local 
source of carriers; G n and Gp represent the rate of carrier generation, and. R n 
and Rp represent the rate of recombination. These factors each depend upon the 
details of the structure of the semiconductor, such as the location and nature of 
imperfections or the distribution of donor and acceptor ions. Information about 
these structural features is scarce at present, mainly because the most successful 
semiconductor-electrode materials are metal-excess semiconductors derived from 
perovskite compounds which have not been well studied. To circumvent this lack of 
detailed information, a series of approximations to the carrier source term are 
described. The most significant of these approximations is the use of a mean carrier 
lifetime which is constant over a variety of different conditions. 

First, consider that carriers are generated in pairs. That is, each carrier is 
produced by the excitation of an electron from the valence band to the conduction 
band. Consequently, a hole is created in the valence band, and the rate of electron 
generation equals that of hole generation. Furthermore, pairs are generated ther- 
mally at a rate given by g Q , and they may be photogenerated by light of intensity 
Iq at a rate of . Thus, we have a net carrier generation rate 


G n ” g 0 + g i Q ” G p (5 * 8) 

The thermal generation rate for electrons and holes in the bulk of the semiconductor, 
where there is no illumination, should equal the corresponding equilibrium recombina- 
tion rates. These rates in turn depend upon the bulk concentrations and the thermal 
equilibrium value of the mean carrier lifetime. The thermal generation rate through- 
out the semiconductor and under all conditions of illumination are given by 


n b /T n " 9 0 ’ V T p 


(5.9) 


Also, assume that recombination of the charge carriers occurs in pairs within the 
semiconductor. No other mechanisms are considered which provide a sink of carriers 
such as trapping and surface recombination. Hie recombination rate is then given by 


R = n/x = p/x = R (5.10) 

n n p p 

Also, in the sequel assume that the carrier lifetimes T n and are constant and 

represent the average values of lifetime for all recombination processes occurring in 
the semiconductor. Combining the results of equations (5.8) through (5.10) yields 
the net carrier sources as 


G - R = 
n n 



(n - n )/t 
b n 


(5.11) 
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and 


G - 
P 



(p - P.)/T 
b p 


(5.12) 


Combining the expressions for the current densities with the continuity equa- 
tions yields the carrier transport equation. This equation is most succinctly given 
in matrix form. To do this, the following terms are defined: 
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The transport equation is then expressed as follows: 


- DV 2 u - • (uf) = F(u) 

ot 


(5.13) 


In this equation, the first two terms describe the diffusion of charge carriers, the 
next term accounts for their drift in the electric field, and the last term describes 
the reaction (generation and recombination) of the carriers. 


5.2 Transport in Semiconductor Electrode 

The general transport equation (5.13) is used to describe the behavior of the 
photoproduced excess charge carriers in a semiconductor electrode. As described in 
section 4, consider a semiconducting cylinder of cross-sectional area A oriented so 
that the liquid junction is perpendicular to the x-axis at x = 0. The concentra- 
tions of free carriers are defined for x > 0 and t > 0 and are given by the 

functions n(x,t) and p(x,t). Within the bulk of the semiconductor, these concen- 
trations have the constant values n fa and p^, respectively. The charge distribu- 
tion in the semiconductor is due to the free carriers and to the distribution of 
ionic charge (assumed to be uniform throughout the semiconductor), so that 

p(x,t) = e [p(x,t) - n(x, t) - (P b - n fo )] (5.14) 
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This charge distribution is the source of an internal electric field E(x,t), as 
determined by the Poisson equation 


-> 

eV 


E = p ( x , 1 ) 


It is this field which is responsible for the drift in equation (5.13). 

The inherent nonlinearity in this system arises from the fact that the field 
which causes the mobile carriers to separate depends upon the local distribution of 
carriers, as shown in equation (5.14). To gain some insight into the effect of this 
nonlinearity, we first examine the sourceless (F = 0) and purely kinetic (D = 0) 
case. In this way information is obtained about the possible stable asymptotic (in 
time) states. Equation (5.13) with F = 0 = D becomes 


■§£-■!« P<u) (5.15) 

As shown subsequently, the concentration for each of the charge species satisfies an 
equation displaying a quadratic nonlinearity which can be approximated by the form 


Z t = a(z 0 ~ Z ^ Z 


(5.16) 


The phase portrait for this ordinary differential equation (shown in fig. 9) displays 
two asymptotic states: z = Zq and z = 0. The asymptotic state z = z Q is the 
only stable one, so any (nonzero) concentration tends to this state. We conjecture 
that this situation pertains to all cases where t > 0, even when sources and diffu- 
sion are present. As a consequence there is a stable asymptotic (steady-state) solu- 
tion to equation (5.15). One way to test this conjecture would be to study the 
transient response of this system for different illuminations. What we shall do 
instead is examine the limiting case of very intense illumination and estimate values 
of the incident intensity which lead to saturation of the photoeffects. 

Under steady-state conditions, the transport equation (5.13) becomes 


V 2 U + D -1 n $ • (|jli) + D _1 F = 0 


(5.17) 


Although unnoticed by many researchers, the system of equations (5.17) and 
(5.14) posseses a first integral. This fact can be used to eliminate the nonlocality 
of the system. The local charge excess distributions are as follows: 


Q + (x) 



dx[p(x) 



(5.18) 
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and 


Q“(x) 



dx[n(x) 


n. 


] 


The local field can be expressed as 

E(x) = E( 0) + -^[Q + (x) + Q~ ( x ) ] 
while the steady-state equation (5*17) becomes 
u” - pE ( x ) B • u' + Au = I(x) 

with 

u = [q + ( x ) , Q (x) ] 

and 



(5.19) 


(5.20) 


(5.21) 
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where 


I (x) = — -^[a e E(0) - Q + (“>) ] + -|p- I Q exp(-ax) 
p L p p 


I n (x) = — t[a n e E(0) - Q (")] - - 5 - I Q exp(-ax) 


n 


2 1 s 

p T 2^ p J 


a = p — P, L 2 
p K e b p 


= rV 1 + O 


" 

n 


a - p i n K L 2 
n r e b n 


Equation (5.17) describes a system with a quadratic nonlinearity and suggests an 
iterative scheme for solving the system. Let u 0 be the solution to 


11 

u Q - Au Q = I(x) 


and from this solution form the drift field 


E(x) = E( 0) + Jl(u ) 


where !( ) is a linear functional. Next, solve equation (5.17) with this drift 
field. In other words, u Q is the first solution in a sequence generated by the 
progression 


u 

n 


u _ + 6u . 
n-1 n-1 


(n — 1, 2, 3, ...) 


where each solution satisfies the linear differential equation 


u n _ Au n = I(x) + Pt E(0) + A ( u n-l) u n-l ] 


(5.22) 
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By analogy with the sourceless kinetic case described by equation (5.15), equa- 
tion (5.22) is expected to possess one stable solution which will be reached after a 
sufficient number of iterations. 


5.3 The Effect of High -In tensity Light 

Examine the first-order solution to the system of equation (5.17) subject to 
boundary conditions at x = 0 of the following form: 


n( 0) = n fe exp(pAV) 


and 


p(°) = P b exp(-pAV) 


where AV is the potential drop across the semiconductor. This solution leads to a 
relationship between the potential AV and the incident light intensity of the form 


^ - exp(-pv) + 
0 



(5.23) 


where 
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0 



b L n 


2 


and 


I 


0, s 


= “Vn 


Equation (5.23) can be rewritten in terms of X = pAV as 


Y~ = \(X f X Q ) = -L + 1 - exp(X) 
0 f s 0 


and it gives the ratio of I Q to Iq q in terms of a single parameter X Q = pv Q . 
Figure 10 depicts this relationship. Ihere are two zeros for \(X,X Q ) : one at 
X = 0 and the other at X = -Xq. Only one of these can correspond to a stable 
asymptotic state, and the well-known behavior of V for low to moderately high 
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intensity suggests X = -X Q . Thus, V Q is identified with the dark bent-band poten- 
tial as 


V = V = V - V 
0 bb F REDOX 


Inverting \(X,X Q ) and taking the stable branch indicated by the solid line in fig- 
ure 9, the bent-band potential is found as a function of light intensity. This curve 
maximizes for 


I 


0 


0 


/ 


s/max 


In X r 


+ 1 


as shown in figure 11. For X » 1 near the saturation point, equation (5.23) can be 
approximated by 


AV « V 




(5.24) 


Thus, an estimate for the intensity of the incident photon flux at which saturation 
occurs is obtained as follows: 


= covD 

0,s bn 


(5.25) 


To make a general estimate of saturation levels, take a = 10 6 m \ n = 1 0 24 m 3 , 
and D = 10” 3 m 2 /s to obtain 


27 2 

I » 10 photons/m -s 
u , s 


This corresponds to an energy flux of 1 0 4 W/cm 2 for monochromatic light at wavelength 
0.55 p. This is a conservative estimate, since the quantum efficiency of the process, 
other charge sinks, and the nonlinearity in the recombination process that would 
require increased flux have been neglected. 


6. SATURATION OF PHOTOELECTROLYSIS OF WATER USING 
A PHOTOELECTROCHEMICAL CELL 

In recent investigations of photoelectrochemical cells with semiconductor elec- 
trodes, there appears to be no evidence of saturation of the energy conversion pro- 
cess (refs. 15 through 18). The photocurrent at a fixed cell voltage increases 
linearly with light intensity up to 380 W/cm 2 (ref. 13). The energy flux level at 
which saturation occurs is an important parameter which is needed in order to under- 
stand the processes governing the operation of the cell. In addition, the saturation 
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flux level is a critical parameter for the design of effective energy conversion 
devices. In the theoretical study of nonlinear effects of charge transport processes 
within the semiconductor electrode , it was found that the light intensity which 
causes saturation varies directly with the charge carrier concentration (eq. (5.25)). 
An experimental search for saturation effects using SrTiO^ crystals was undertaken. 
These crystals were slightly reduced to form a metal excess semiconductor of high 
resistivity (low carrier concentration) . The SrTiO^ electrodes were illuminated with 
the 351 -nm and 364-nm radiation from an argon-ion laser with a maximum intensity of 
60 mW/cm^. A nonlinear variation of the photocurrent with light intensity and the 
evidence for complete saturation of the photoprocesses in SrTiO^ a light intensity 
near predicted values were observed. This experiment is discussed in the following 
paragraphs • 

Single-crystal wafers of SrTiO^ were cut from a boule and polished. The grit 
size of the final polish was 0.05 jam. The crystals were reduced in a hydrogen atmo- 
sphere at about 700° C for less than 1 minute. The conductivity was estimated by 
measuring the current as a function of applied voltage for a fixed electrode separa- 
tion using a curve tracer. An attempt to determine the resistivity of the samples by 
the four-point probe technique (ref. 19) proved unsuccessful because of high sample 
resistivity. Ohmic contacts were made using an indium alloy rubbed on one side of 
the crystals and a copper wire attached to this contact with a silver paste. The 
sample was mounted on a glass tube with a silicone-based adhesive. The electrolytic 
cell was a stainless -steel vessel with a quartz window and consisted of the SrTi0 3 
photoanode , a platinum counter electrode, a saturated calomel reference electrode 
(SCE), and a buffered solution of NaOH adjusted to pH-13. During the experiment, 
hydrogen gas was bubbled across the platinum electrode and, using a potentiostat, the 
photoanode was maintained at zero volts with respect to the SCE. The level of uv 
radiation from the argon-ion laser was controlled using neutral density filters cali- 
brated for uv radiation. The laser output was continuously monitored by reflecting 
part of the laser beam into a power meter. 

Thirty-two crystal samples were prepared as described above and used to form the 
photoanode of an electrolytic cell. Characteristic current-voltage data were taken 
at fixed light intensity values to insure that each cell was operating properly. 

Then, while the cell remained at a fixed bias of zero volts relative to the SCE, the 
intensity of the incident radiation was varied, and the photocurrent was measured for 
each light intensity. The results of these measurements separated the crystal sam- 
ples into three categories, distinguished generally by crystal resistivity. Three of 
the samples had a very high resistivity (on the order of 1 0 6 ohm-cm), and each of 
these exhibited a saturation in photocurrent with light intensity. A typical satura- 
tion curve for a crystal in this category is shown in figure 12. Several samples had 
resistivities varying from 1 0 4 to 10 3 ohm-cm. For these samples the photocurrent 
versus intensity displayed some curvature, suggesting saturation but not achieving it 
at the intensities available in this experiment. The remaining category of crystals 
had resistivities of 10^ ohm-cm, and for these samples the photocurrent varied 
linearly with light intensity within the intensity range of this experiment. 

The model of excess charge transport processes within the semiconductor devel- 
oped in the preceding section predicts a saturation for the photon flux given by the 
expression 
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where is the bulk free carrier concentration, D n is the diffusion constant, 

and a is the photoabsorption coefficient. Since the semiconductor resistivity p 
is related to the concentration n fa and mobility p = pD (where p = e/kT) of the 
carrier by p = (en^li^) ' equation (6.1) can be written as 
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0,s 
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e PP 


( 6 . 2 ) 


Values for the absorption coefficient of SrTiO^ for radiation in the uv region are 
not readily available, but an average value of 4.4 x 1 0 4 cm” 1 is quoted by Capizzi 
and Frova (ref. 20). For a crystal with a resistivity of 1 0 6 ohm-cm and an inci- 
dent photon energy equal to that of the argon laser used in this experiment, equa- 
tion (6.2) gives a saturation flux of 2 mW/cm 2 . This value agrees with the order-of- 
magnitude estimate for the experimental results presented here. A fully reduced 
crystal of SrTiO^ has a resistivity about 10 6 times lower than the samples used here, 
so photosaturation at an energy flux of 2 kW/cm 2 would be expected. This high value 
of radiant energy flux is consistent with the failure to observe saturation in pre- 
vious experiments • 


7. THERMAL PHOTO ELECTROCHEMICAL CELL 

The solar conversion efficiency of a stable semiconductor-based photoelectroly- 
sis cell is on the order of 1 percent. This is a result of the fact that only the 
metal-oxide semiconductors are stable to the electrolysis process, and they are pho- 
toactive for photon energies around 3 eV or for uv radiation. Decreasing the bandgap 
of the stable electrodes should increase the efficiency, since more of the solar flux 
would be utilized. The bandgaps of most semiconductors decrease with increasing 
temperature, so an experimental thermal photoelectrochemical cell was assembled to 
assess the effect of temperature on conversion efficiency to broadband radiation. 

Two temperature effects were anticipated - thermal bandgap narrowing of the photo- 
electrode and decrease of the potential required for the electrolysis of water 
(reversible electrolysis potential). Both of these effects should result in an 
increase in the solar conversion efficiency. Figure 13 shows the effect temperature 
has on thermoneutral potential and the reversible electrolysis potential for water. 

It can be seen that the reversible electrolysis potential decreases from 1.23 V at 
standard conditions as the cell temperature increases. This decrease is due to the 
increased contribution to the electrolysis reaction made by the thermal energy avail- 
able in the aqueous solution. 

A schematic of the experimental setup is shown in figure 14. The thermal photo- 
electrochemical cell was a stainless-steel vessel capable of providing temperatures 
to 200°C. A quartz window mounted on the side of the stainless-steel vessel allowed 
the photoelectrode to be illuminated by a 150-W xenon lamp. A platinum screen was 
used as a counter electrode, and mercury-mercuric oxide was used as the reference 
electrode. Argon was used to pressurize the vessel to prevent the water from boiling 
at high temperatures. The 0.1 -M NaOH electrolyte and the electrodes were contained 
in a quartz beaker within the stainless-steel pressure vessel. The solution was 
stirred using a Teflon-coated magnetic stirring bar, and the photocurrent-voltage 
characteristics of the cell were monitored with a potentiostat and x-y recorder. The 
photoelectrodes tested were Ti0 2 and SrTiO^ single crystals reduced in a hydrogen 
atmosphere for 2 hours at 800°C. The temperature-induced changes in the photo- 
response characteristics of an SrTiO^ based cell are shown in figure 15. Two thermal 
effects are evident. First, the magnitude of the onset potential (the most negative 
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potential which produces current) decreases as the cell temperature increases. This 
decrease is deleterious to the conversion efficiency, because it is necessary for the 
onset potential to be negative with respect to the redox potential of hydrogen. 
Second, the photocurrent increases with increasing temperature for cell voltages 
above about -0.2 V referenced to Hg/HgO. This effect favors the increase in the 
conversion efficiency. The conversion efficiency r\ can be obtained from the 
relation 


I . (l .23 - V ) 

= ph a PP ; 

^ Photon flux 


(7.1 ) 


where I * is the photocurrent, v app is the externally applied voltage, and 1.23 
is the reversible electrolysis potential for water. A relative efficiency can be 
determined by using the efficiency at room temperature as a reference. A plot of the 
relative efficiency as a function of cell temperature for an SrTiO^ based photoelec- 
trochemical cell is shown in figure 16. Hie relative efficiency peaks at about 100°C 
and shows a 50-percent increase in efficiency. Hie decrease in conversion efficiency 
above 100°C is caused by the thermally induced decrease in the onset potential for 
this electrode. The relative efficiency for the SrTiO^ samples peaked at tempera- 
tures ranging from 75°C to 100°C. Hie thermal response of TiC> 2 based cells was simi- 
lar to the SrTiO^ based cells, that is, onset potential decreased with increasing 
temperature, photocurrent increased with temperature, and the relative efficiency 
peaked at temperatures ranging from 75°C to 100°C. The rate at which the onset 
potential decreased with increasing temperature for the TiC> 2 samples was smaller than 
for the SrTiO^ crystals. This effect was investigated and found to be related to the 
point of zero zeta potential (pzzp) of the material. A discussion of the pzzp 
results is included in the next section. 

The relative photocurrents for both Ti0 2 and SrTiO^ at 1.5 V referenced to 
Hg/HgO as a function of temperature are shown in figure 17. The thermal increases in 
the relative photocurrent for both the SrTiO^ and TiC> 2 correspond to approximately 
0.25 percent/K. The thermal enhancement of the photocurrent observed here is attrib- 
uted to changes in the photoelectronic properties of the semiconductor anodes. The 
absorption coefficient of the metal oxides is known to be temperature dependent 
(refs. 20 through 22) and is considered to be the principal cause of the photothermal 
enhancement observed. 

The spectral dependence of the absorption coefficient a for the two semicon- 
ductor materials used can be divided into two energy regions. One energy region is 
below the bandgap and is governed by the Urbach rule (refs. 20, 21, and 23), and the 
second region of energies is greater than the bandgap energy. 

Figure 18 illustrates the two absorption regions as well as the spectral irradi- 
ance of the xenon lamp, which is nearly constant in the region between the bandgap of 
SrTiO^ (400 nm) and 420 nm. The Urbach effects for the two cell-temperature extremes 
(300 K and 450 K) are shown here and are calculated from 
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where a Q = 4,4 x 10 \m (ref. 18), T is the absolute temperature, E^ is the 
bandgap, k is the Boltzmann constant, h is the Planck constant, c is the speed 
of light, and X is the wavelength. "The temperature dependeng^ of the current dens- 
ity produced by absorption in the Urbach region is eFI o(h”/ w ^ ere is the 

light intensity at the bandgap and F is a constant of order unity. Even though 
there is a significant change in the absorption coefficient in the Urbach region due 
to temperature changes, this effect contributes only approximately 1 percent of the 
total thermal-induced enhancements observed. 

Assuming the absorption coefficient is independent of changes in bandgap, the 
change in the relative photocurrent with temperature is approximately 1 percent/K. 

To arrive at this figure, the spectral output of the light source and the rate of 
change of the bandgap with temperature for SrTiO^ (-0.95 x 10 -4 eV/K) were used 
(ref. 24). The results for TiO^ are expected to be similar. This result is a factor 
of 4 higher than observed experimentally, but no consideration has been given to the 
effect of temperature on loss mechanisms such as recombination. 

In conclusion, experiments indicate that photocurrent efficiency of broadband 
radiation for SrTiO^ and Ti0 2 based photoelectrochemical cells can be enhanced and 
that this increase can be attributed to the thermal decrease in bandgap. This tech- 
nique may also be beneficial in enhancing the solar conversion efficiency for the 
stabilized liquid junction solar cells. 


8. ONSET POTENTIAL AND POINT OF ZERO ZETA POTENTIAL 

In the previous section it was noted in the discussion of figure 15 that the 
onset potential for both Ti0 2 and SrTiO^ samples decreased as the cell temperature 
increased, and that the rate observed for Ti0 2 was smaller than that for SrTiO^. The 
thermal change in the onset potential for single crystals of a-Fe 2 0 2 and a polycrys- 
talline WO^ sample were also investigated. 

The onset potential was taken from the photocurrent voltage recordings using a 
conventional three-electrode system. The experimental apparatus is shown schemat- 
ically in figure 19. Ihe light from a 1 50-W xenon lamp was chopped and the photo- 
current recorded using a lock-in amplifier. The cell was held inside a stirred and 
heated water bath. A reversible hydrogen electrode (RHE) was used as a reference 
electrode for most of the experiments. In some experiments conducted in basic solu- 
tions, an Hg/HgO electrode was used as a reference electrode for comparison. Solu- 
tions of NaOH (pH-13) and H 2 SC> 4 (pH-2) were used. 

The effect of temperature on the onset potential for these materials at two 
different solution pH values is shown in figure 20. A linear change of the onset 
potential with respect to the RHE is observed in all the cases. The slope of the 
change is different for each electrode material; the highest slopes were observed for 
Fe 2°3 anc ^ SrT i-° 3 / and the lowest for WO^. These measurements were repeated with 
solutions of different pH. From figure 20 it can be seen that the slope of the onset 
potential versus temperature is about the same for each electrode with solutions of 
either pH-13 or pH-2. 

The relationship between the onset potential (taken to be the flatband poten- 
tial V FB for the electrode), the redox level for hydrogen V Q , the potential drop 
across the Helmholtz layer V H , the electron affinity EA, and the position of the 
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Fermi level with, respect to the conduction band edge is shown in figure 21 and 

is given (refs. 25 through 28) as 


EA = eV Q + eV H + eV pB + A pc (8.1) 

Changes in temperature can effect A , V H , and Vq. It can be shown (ref. 29) 
that the change in the onset potential with temperature using an RHE as the reference 
electrode is 


AV 

FB = 2 . 3R H 
AT F P pzzp 



( 8 . 2 ) 


where R is the gas constant, F is the Faraday constant, n is the free carrier 
density in the semiconductor, N c is the density of states in the conduction band, 
and 


PH 


pzzp 
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which is the pH at the point of zero zeta potential for the electrode. Equa- 
tion (8.2) shows that the onset potential of an oxide semiconductor electrode shifts 
linearly with increasing temperature to a more positive potential with a slope which 
is independent of the pH of the solution and which is determined by the magnitude of 
the pH . The more basic the electrode surface, the higher the and the 

PZZp pzi p 

result is a larger change in the onset potential with cell temperature. The P H p ZZ p 
for each electrode was determined using equation (8.2) and the slopes of the curves 
in figure 20 for pH-2 and pH-13 solutions. These values are given in table I. Also 
given in table I are the carrier concentrations and density of states for each of the 
semiconductor materials used to determine the magnitude of the second term on the 
right-hand side of equation (8.2) • The average value for the P H p ZZ p of each elec- 
trode is shown in figure 22. It should be noted that the tungsten trioxide sample is 
degenerately doped, and therefore the second term on the right-hand side of equation 
(8.2) is zero. Also shown in table I and figure 22 are P H p ZZ p values for these 
semiconductor materials taken from the literature. The curve arawn in figure 22 is 
taken from reference 26. There is good agreement between the values obtained in this 
work and the values taken from the literature. Since equation (8.2) properly 
describes the observed behavior of the onset potential of the oxide semiconductors 
reported here, the changes in the onset potential with temperature are attributed to 
thermal effects on the potential across the Helmholtz layer. 

These results lead to two important considerations for choosing semiconducting 
materials to be used as photoelectrodes in a photoelectrochemical cell. First, the 
onset potential of the photoelectrode is dependent on the EA of the electrode mate- 
rial and on the potential across the Helmholtz layer. The onset potential is more 
negative the smaller the EA and the more acidic the electrode surface (smaller 
pH ). While the EA of the electrode material cannot in principle be changed, 
the 2 Z pH can be changed (ref. 30). It was suggested (ref. 30) that by binding 

more acidic groups to the surface of TiC^, for example, the onset potential could be 
shifted to more negative values. The method described here can be used to evaluate 
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in a relatively easy way any change in the acidity of the electrode surface following 

treatment. (The common methods of evaluating the pH of the materials in the 

• pz z p 

powder form are not applicable here.) 

Second, as discussed in the preceding section, there is the possibility of 
increasing the solar-to-chemical energy conversion efficiency of a photoelectrochemi- 
cal cell by increasing the cell temperature (refs. 28 and 29). For SrTiO^ and TiO^, 
the photocurrent increases by increasing the cell temperature, and the overall effi- 
ciency begins to increase. At cell temperatures higher than about 80°C, the effi- 
ciency begins to decrease because the onset potential for these materials has shifted 
to more positive values. 

These experiments suggest that by decreasing the value of the P H DZZ n ( e *9*/ by 
binding acidic groups to the surface), the amount of the thermally induced shift in 
the onset potential of the photoelectrodes may be reduced. This reduction would 
improve the efficiency of the photoelectrochemical cell by increasing the 
temperature. 


9. PHOTOPRODUCTON OF HALOGENS 


To obtain the maximum conversion efficiency for a given semiconductor, and 
therefore for a fixed bandgap, it is necessary to find a redox couple whose elec- 
trolysis potential is slightly less than the bandgap energy and whose redox reactions 
lie within the bandgap. For example, the maximum efficiency attainable (as shown in 
section 7) for the electrolysis of water (electrolytic potential of 1 .23 V and 
enthalpy of 1.47 V) by TiC> 2 with a bandgap of 3 eV is 


_ 1 .47 eV 
\ax 3.0 eV 


0.49 


(9.1) 


The conversion efficiency can be increased, for example, by finding an electrolytic 
reaction which is driven at about 3 V. In this situation the conversion efficiency 
would be 1 • These efficiencies presume that the incoming radiation is monochromatic 
and has an energy of 3 eV, which is the bandgap of TiC^. 

An example of one such reaction is the electrolysis of HF. In acid solutions 
the electrolytic potential is 3 V, and the reactions at the electrodes produce hydro- 
gen and fluorine as 


2HF -> H + F^ 
2 2 


(9.2) 


Experiments were attempted using SrTiO^ and Ti0 2 single crystals in aqueous solutions 
of 0.1 M HF. However, as expected from the corrosive nature of HF, the HF chemically 
attacked the semiconductor electrodes. Attention was turned to the other halogens: 
chlorine, bromine, and iodine. The redox levels for these halogens and water in 
relation to the bandgap of TiC> 2 at two acid pH values are shown in figure 23. The 
photoelectrolysis of these compounds has been analyzed (refs. 31 and 32) for the case 
of single-crystal Ti0 2 photoelectrodes in the standard photoelectrochemical cell. 

The results of these analyses indicated that chloride, bromide, and iodide ions were 
oxidized at an illuminated Ti0 2 electrode. 
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Some work has been done using n-type semiconductor powders (refs. 33 through 37) 
to photocatalyze oxidaton reactions which otherwise would occur at very low rates. 

The platinized semiconductor powder can be considered to provide many short- 
circuited, platinum, semiconductor electrodes on which photoassisted reactions can 
take place, provided the reduction process occurs at a more positive potential than 
the flatband potential of the semiconductor. To date, most of the reactions at illu- 
minated TiC^ or Pt-TiO^ powders are, by their nature, photocatalytic (refs. 33 
through 39). The possibility of carrying out photosynthetic reactions is of greater 
interest (refs. 39 and 40). The experimental results summarized herein deal with the 
producton of the halogens at TiO^ powders suspended in aqueous solutions of their 
respective ions. The reactions investigated were the photoassisted oxidation of the 
halogens and the reduction of oxygen at TiC> 2 particles. The process is shown sche- 
matically in figure 24 for chlorine. The two half -reactions assumed are: 

2X~ + 2h + * X 2 (9.3) 

at the illuminated Ti0 2 site, and 

J 0 2 + 2H + + 2e~ -► H 2 0 (9.4) 

at the platinized site. These half -reactions give the overall reaction 

2X " + \ °2 + 2H+ * X 2 + H 2 ° (9#5) 

where H + and e“ are the photoproduced hole-electron pair and X is the appro- 
priate halogen. The experimental details can be found in reference 41 . 

The energetics of the processes described by reactions (9.3) through (9.5) for 
the three halogens are shown in figure 23. The figure shows that the oxidation of 
1“ to by oxygen is a spontaneous reaction. Adding the Pt-Ti0 2 photocatalyst, 

however, was found to enhance the oxidation rate, as shown in table II. More strik- 
ing is the photo-oxidation of Br” and Cl” at the Pt-Ti0 2 powder surface. The oxida- 
tion of these halides by oxygen is nonspontaneous at pH - 5.5 with AG = 8 kcal/mole 
for Br oxidation and AG = 21 kcal/mole for Cl* oxidation. However, these halides 
were produced in the presence of the platinized titanium dioxide particles and light, 
and the results are given in table II. These are clear examples of a chemical reac- 
tion taking place at the illuminated platinum-semiconductor powder surface in a non- 
spontaneous (photosynthetic) direction and giving rise to useful chemicals. 

The fact that photosynthesis is clearly achieved in this system indicates that 
charge separation occurs in the semiconductor powder (without external biasing) by 
introducing the platinum. With pure powder, much smaller rates were found, indi- 
cating charge separation was not effective. (See table II.) The effect of platini- 
zation is associated with the fact that the overvoltage for oxygen reduction at the 
particle is lowered (refs. 38 and 42). 

It was found that the rate of production of the halogens is in the order 
*2 > Br 2 > C ^2* is because of the difference in the ability of the halide 
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anions to compete with water for the photoproduced holes in the semiconductor 
(refs, 31 and 32), and also because of the difference in the ability of oxygen to 
compete with the halogen for reduction on the platinum site. 

The rate of the halogen production was enhanced by lowering the solution pH. 

(See table II.) This can be understood from reactions (9.3) through (9.5) and fig- 
ure 23: (a) the oxidation potential of water is more positive at lower pH; there- 

fore, halide oxidation becomes more competitive; and (b) the reduction potential of 
oxygen is more positive at lower pH; therefore, oxygen becomes more competitive with 
the halogen molecules for reduction at the platinum sites. In fact, at the lower pH, 
the oxidation of Br~ becomes photocatalytic as can be seen from figure 23. 

For the case of the photosynthetic reactions demonstrated, net conversion of 
solar energy to chemical energy is also achieved. The efficiency for this conversion 
can be calculated using the equations 


Energy stored as - Energy from power supply 
Light energy 


(9.6) 


ifv . - V It 
v prod ext J 

Light energy 


(9.7) 


where t is the illumination time, I is the rate of X 2 production expressed in 
units of current, V rod is the difference in the redox potentials of the two half- 
reactions of the cel?y°and v ext is th e voltage supplied from an external source. 

The photosynthetic reaction is carried out without supplying electrical power from an 
external source; thus, v ext = 0- Taking into consideration only the light absorbed 
by the TiC> 2 , the calculated optical power efficiency for the ideal case of Cl 2 pro- 
duction is 3 percent. This is compared with the 1- to 2-percent solar-to-electrical 
energy conversion efficiency at TiO^ electrodes (refs. 43 and 44). In the acidic 
solution where Cl^ is allowed to accumulate, a power efficiency of 0.03 percent is 
calculated for the conversion of solar-to-chemical energy. The practical values for 
the efficiency would be between these values. If the whole spectrum of available 
light energy is taken into consideration, the efficiencies calculated are lower. 

The halogens produced in this way can be used not only as a source of chemical 
energy in space and on Earth, but also by the chemical industry, which uses these 
chemicals in large quantities. At the present time, these chemicals are being pro- 
duced primarily from electrolytic cells using electrical energy to drive the desired 
reaction. This technique only requires a source of uv radiation which can be sup- 
plied by the Sun. 

Figure 25 is a schematic of a photoelectrochemical cell as well as a magnified 
view of the chemical processes taking place on the powder. Solar ultraviolet light 
absorbed by the platinized Ti0 2 powder oxidizes the halide ion to the halogen mole- 
cule (e.g., chloride ion from a brine solution to chlorine gas) and reduces oxygen 
which is supplied by bubbling air through the reaction vessel. The chemical reaction 
taking place in the presence of the powder and uv light is shown in reaction (9.5). 
The halogen gas is collected above the brine solution and may be separated from the 
air by membranes. The TiO^ powder is inexpensive because it is widely used in paint, 
and a natural source of brine on Earth could be seawater. The technique of getting 
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chlorine from seawater is an energy-efficient process capable of producing a product 
widely used in the chemical industry. 


10. LAYERED TRANSITION METAL THIOPHOSPHATES 

The search for stable, small-bandgap semiconductor materials for use as elec- 
trodes in photoelectrochemical cells has led to the metal chalcogenides , MX 2 , such as 
MoS 2 , and WSe 2 (refs. 45 through 51). These materials have a layered crystal struc- 
ture. In such structures, the chalcogenide bonds within the layer are saturated, and 
the individual layers are bound to each other by weak van der Waals forces. 
Electron-hole generation appears to occur because of d-d phototransitions. These 
transitions do not affect the chemical bonds in the crystal, so the crystal resists 
photodecomposition. 

A new class of layered compounds was recently analyzed for use as photoelec- 
trodes in photoelectrochemical cells (ref. 52). These materials are of the type MPX^ 
where M is a transition metal, P is phosphorous, and X is an element of the sulfur 
series. These materials have been called transition metal thiohypophosphates 
(refs . 53 and 54) . 

Single crystals of the transition metal thiophosphate crystals were grown by 
iodine vapor transport in quartz ampules 10 cm in length and 1.5 cm in diameter 
(refs. 53 and 54). A 4-g charge consisting of a stoichiometric mixture of the ele- 
ments, to which 150 mg of iodine was added, was introduced into one end of the 

ampule, and the ampule was sealed in a nitrogen atmosphere. A tube oven with a hori- 
zontal temperature gradient given in table III was used for the growth process. 

Growth conditions were maintained for about 100 hours, resulting in crystals up to 
about 10 mm x 10 mm x0.1 mm. 

The photoelectrochemical response of SnPS^, NiPS^, and FePS^ used as electrodes 
in an aqueous solution containing sulfuric acid at pH-2 is shown in figure 26. The 
curves labeled (a) are the dark currents, and the curves labeled (b) indicate the 
photocurrent when the sample was irradiated with the full intensity of the 150-W 

xenon lamp. Ihe three materials shown have a very low dark current and a clear pho- 

toresponse. The SnPS 3 sample has a photoanoidic response typical of n-type semicon- 
ductor compounds. Single-crystal samples of SnPS^ grown under apparently similar 
conditions showed photocathodic response as well as a response expected of an intrin- 
sic semiconductor electrode. This suggests that the growth conditions are important 
in determining the doping characteristics of the SnPS^. Also, two types of SnPS^ 
crystals were consistently produced. One type was distinctly layered with, only a few 
layers per sample, and the other type was more granular and layered. No effort was 
made to correlate the growth conditions with the resulting doping characteristics of 
the SnPS 3 crystals. The FePS^ and NiPS^ crystals behave as p-type semiconductor 
electrodes. Crystals of FePS^ and NiPS^ from different growths all behaved as p-type 
materials • 

Crystals of SnPSe 3 and MnPS^ showed n-type response. These two materials were 
very unstable in both acid and base solutions. The photoresponse of BiPS 4 crystals 
was typical of an intrinsic photoelectrode. This material was also unstable in acid 
or base electrolytes, and as a result, no bandgap data were obtained. In addition to 
the above materials, CoPS 2 , CrPS^, FePSe^, and GaPS 2 crystals were grown and gave no 
photoresponse and were very unstable in acid solutions. 
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The onset potential in volts for the photocurrent of the SnPS^ samples was 0 
referenced to SCE at pH-2. This onset potential is about 850 mV more negative than 
the reversible potential for oxygen evolution, which indicates a net gain in energy 
conversion. The onset potential for FePS^ is about 400 mV more cathodic than the 
reversible potential for hydrogen evolution. This indicates a net gain in energy 
conversion. The position of the band edges of SnPS^ and iodine redox potentials are 
shown in figure 27. 

Investigations of the stability of these compounds under illumination were made 
in both acid and base electrolyte. The crystals were very unstable in a high pH 
environment, and after a short period of time, the photocurrent decayed to zero. 
Inspection of the crystals indicated that delamination of the layers occurred as well 
as photocorrosion. All the selenide compounds photocorroded in acid solutions also. 
The sulfur-based compounds were found to be stable for periods up to an hour in 
strongly acid solutions. The surface of these materials exposed to the electrolyte 
showed no signs of degradation when inspected under a light microscope. However, 
even the sulfur-based crystals were very unstable in acid solutions if the edges of 
the layers were exposed to the electrolyte. 

The photocurrent response as a function of photon energy was determined for some 
of the layered materials. The photocurrent response for n-type SnPS^ is shown in 
figure 28. Assuming that a Schottky barrier exists at the electrode-electrolyte 
interface, then, near the band edge, the photocurrent I , will have the following 
dependence on the photon energy E ph : 



where A is a constant for a fixed applied potential, E G is the bandgap energy, 
and n depends on the bandgap optical transitions. That is, n = 1 for a direct 
transition and n = 4 for an indirect transition. (See refs. 55 and 56.) A plot 
of (E Yi^ph^ ^ as a f unct i- on °f , E ph. for each of the materials examined showed a 
linear ^dependence, which indicated an indirect bandgap. An example of this depend- 
ence is shown in figures 29 and 30 for n-type SnPS^ and p-type FePS^. The intercept 
was the bandgap of the material. The bandgaps determined by this method for these 
crystals are listed in table III. 

Experiments to determine the effect of solution pH on the onset potential showed 
a linear dependence up to PH-6. Figure 31 gives the data for SnPS^. The curve has a 
slope of 59 mV/pH and indicates approximately Nernstian behavior of the onset poten- 
tial for these layered materials. 

The spectral response of the SnPS^ photoelectrodes was investigated as a func- 
tion of cell temperature in a pH-2 solution. The layered sample exhibited a broad 
response centered at 3 eV with a full width at half maximum of 1 eV as shown in fig- 
ure 28. The granular samples have a much narrower response centered at 2.3 eV with a 
full width at half maximum of 0.15 eV as shown in figure 32. The effect of cell 
temperature on the spectral response of a granular SnPS^ sample is also shown in 
figure 32. There is a dramatic shift of the response toward the red end of the spec- 
trum with increasing temperature. A plot of the peak of the spectrum as a function 
of cell temperature is shown in figure 33. The shift of the peak is at a rate of 
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2.4 meV/K. This thermal shift, however, is not due to the change in the bandgap, 
since no significant change in the bandgap («2.1 eV) was observed. 

Hie broadband photocurrent efficiencies of the SnPS^ samples were on the order 
of 0.01 percent. The monochromatic quantum efficiency for light corresponding to the 
peak in the spectral responses of the two types of SnPS^ crystals (2.3 eV for the 
granular SnPS^ and 3.0 eV for the layered SnPS^) was 8 percent and 40 percent, 
respectively • 

Experiments were conducted using the layered crystals as a photosensitive elec- 
trode in a liquid-junction photovoltaic cell using a redox couple made up of iodide 
and tri -iodide in 1 M sulfuric acid. The materials previously found to be unstable 
were also unstable in the iodide-acid solution. The experimental photocurrent- 
voltage curve for n-type SnPS^ in a solution of 1 M H 2 SC> 4 , 0.025 M I 3 , and 1 M I* is 
shown in figure 34. Light from the 150-W xenon lamp was focused on the electrode. 

The open-circuit voltage was about 0.31 V, the fill factor was 0.38, and there was an 
efficiency of 4.5 x 10~ 2 percent. Other crystals showed low fill factors as well as 
low conversion efficiencies. Hiis may be a result of surface recombination losses, 
high absorption of the incident radiation by the iodide solution, and high resistance 
exhibited by the layered materials. 

In conclusion, a new class of layered transition metal thiophosphates of the 
type MPX3 were investigated for use as photoelectrodes in photoelectrochemical cells. 
Only SnPS^, NiPS^, and FePS^ crystals were stable, and these were stable only in acid 
solutions. Crystals of SnPS^ exhibited intrinsic and n- and p-type photocurrent 
characteristics. The iron-based and nickel-based sulfur compounds showed only p-type 
photocurrent-voltage characteristics. These materials had an indirect bandgap of 
about 2.1 eV and an onset potential near zero referenced to the SCE at pH-2. The 
onset potential for the stable materials was Nernstian between pH-0 and pH-7. These 
materials showed low conversion efficiencies when used as photoelectrodes in liquid- 
junction photovoltaic cells with iodide and tri-iodide as a redox couple. Based on 
these results, further research characterizing the growth conditions resulting in n- 
and p-type crystals, particularly for the sulfur-based crystals, as well as more 
extensive investigations with other redox couples will determine the utility of this 
class of compounds as photoelectrodes in photoelectrochemical cells. 


1 1 . HIGH-EFFICIENCY PHOTOELECTROCHEMICAL SYSTEM 

Laser energy conversion efficiencies of 20 percent for the photoelectrolysis of 
water to gaseous hydrogen and oxygen in a SrTiCU based cell with a small applied 
voltage have been demonstrated by Wrighton et al. (ref. 16). The calculation of the 
conversion efficiency must account for this applied potential as 



( 1 .48 - V 


app 


Light flux 


) 




where 1.48 is the reversible potential corresponding to the enthalpy for the reaction 


2H 2 + °2 ♦ 2 H 2 ° 


( 11 . 2 ) 
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Elimination of the applied potential from equation (11.1) increases the conversion 
efficiency, and, for the case of a light flux from a laser, the efficiency is limited 
by the quantum efficiency of the particular photoelectrode. Approaching a 40-percent 
conversion efficiency for this SrTiO^ system is feasible. 

One scheme which successfully demonstrates that high laser conversion efficien- 
cies can be achieved is shown in figure 35. Laser energy incident on a photoelec- 
trode must pass through a window in the photoelectrochemical cell. A small portion 
of this laser energy is reflected from the window and would normally be a loss in the 
conversion process. However, placing a photovoltaic device, such as a Schottky bar- 
rier solar cell, in the reflected beam as shown in the figure can provide the 
required external bias necessary to achieve saturation and maximum efficiency. 

A system, configured as shown schematically in figure 35, consisting of a single 
crystal of a-Fe 2 0^ as a photoanode in the photoelectrolysis cell and three silicon 
solar cells was used to demonstrate the principle of high laser-to-chemical energy 
conversion efficiency. The open-circuit voltage of the silicon solar cells with 
laser radiation having a photon energy of 2.4 eV was about 0.3 V. Figure 36 is the 
experimentally obtained photocurrent response of the iron-oxide crystal as a function 
of voltage. The response is shown as follows: (1) Unilluminated; (2) Illuminated 

with laser radiation at E , =2.4 eV; (3) Illuminated using as a bias the portion of 

the laser radiation reflected from the cell window and surface of the iron-oxide 
crystal. The onset of the photocurrent shifts approximately 1 V with respect to the 
hydrogen half -cell potential. The principle of the technique is, therefore, verified 
by this result. Not only is the loss inherent in the reflection from the surfaces 
recovered, but also the photovoltaic cells utilize this energy to provide the bias 
needed to maximize the conversion efficiency by eliminating the need for externally 
derived power. 

The current efficiency of a photoelectrochemical cell is also a function of 
externally applied voltage. This is particularly true of indirect bandgap materials 
such as ct-Fe 2 0 3 shown in figure 37. Peak current efficiency occurs for an external 
bias of 2.0 V at a laser photon energy of 3.5 eV. The bottom curve in figure 37 
shows the maximum laser conversion efficiency as a function of laser photon energy in 
eV. The efficiency for converting laser radiation to hydrogen and oxygen peaks at 
3.5 eV and is 15 percent if the 2.0-V bias is not accounted for. 

The top curve of figure 38 gives the log of the quantum efficiency as a function 
of photon energy for single-crystal SrTiO^. It should be noted that the quantum 
efficiency is unity for photon energies near 3.8 eV. Since the technique as demon- 
strated eliminates the bias previously required to achieve saturation current, laser 
conversion efficiency for this material and any other material is now determined by 
its spectral response. In the bottom curve of figure 38, the laser-to-chemical con- 
version efficiency is plotted as a function of the incident laser photon energy. The 
efficiency reaches 38 percent for laser photon energies of 3.8 eV. 

Single-crystal photoelectrodes give the highest conversion efficiencies. Fig- 
ure 39 shows that a polycrystalline material (WO^) can also have reasonably high 
conversion efficiencies. This material is made by starting with metallic tungsten 
and then heating the metal in air at about 500°C until the oxide forms. The process 
is simple and therefore inexpensive. The laser-to-chemical conversion efficiency for 
this material reaches 6 percent at a laser photon energy of 2.8 eV. 
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Table IV shows the semiconductor materials that are stable during the photo- 
electrolysis process. Also included are their respective flatband potentials with 
respect to a reversible hydrogen electrode and the bias needed to give the maximum 
conversion efficiency. The last two columns give the energy of the laser photons and 
the maximum conversion efficiency. In every case, a bias of at least 1.5 V is 
required to achieve the highest efficiencies, and the energy of the required laser 
photons is in the ultraviolet. A Schottky barrier zinc-sulfide photovoltaic cell has 
its peak response in the uv and is capable of supplying the bias required for maximum 
efficiency. (The efficiencies cited herein are based on the enthalpy or energy 
evolved in the chemical recombination of hydrogen and oxygen to form water as would 
be the case, for example, for the combustion process taking place in a rocket 
engine • ) 

The highest efficiency, 38 percent, is obtained by an SrTiO^ electrode, and the 
next highest efficiency (33 percent) is obtained by a TiC> 2 electrode. A conversion 
efficiency of 6 percent can be achieved for a polycrystalline electrode of thermally 
grown tungsten trioxide. 

Unsuccessful efforts were made to construct a ZnS photovoltaic cell to demon- 
strate the high conversion efficiencies. Therefore, the premise for high laser con- 
version efficiency with stable photoelectrodes is not fully demonstrated. 


12. CONCLUDING REMARKS 

An experimental and theoretical investigation of the semiconductor electrode has 
been carried out, and particular emphasis has been placed on the possibility of its 
use as part of a space power system. Using a semiclassical analysis, a model of the 
physical and chemical processes that occur when the electrode is illuminated has been 
developed. Ihis model describes the static space-charge distribution within the 
semiconductor and the transport of photogenerated excess charge to the electrode- 
electrolyte interface, where it is transformed by a chemical reaction. Using an 
a PP rox iroation to the charge transport equations, an estimate was made of the inten- 
sity of incident illumination at which the photogenerated carrier transport would 
saturate. Using SrTiO^ electrodes, the saturation of photoeffects was observed with 
an argon-ion laser. 

Experimentation on the temperature dependence of the broadband photocurrent for 
SrTiO^ and TiO^ photoelectrochemical cells demonstrates an increased efficiency with 
increasing temperature. This effect, which is attributed to the fact that the band- 
gap narrows with temperature, may prove useful in enhancing the solar conversion 
efficiency for liquid-junction solar cells. 

The effect of temperature on the onset potential for the photocurrent was inves- 
tigated for the metal-oxide electrodes Fe 2 ° 3 ' TiC^, SrTiO^, and WC> 3 . In aqueous 
electrolytes, the onset potential decreased linearly with increasing cell tempera- 
ture. An analysis indicates that the slope of this curve is proportional to the 
point of zero zeta potential for the particular electrode material. The point of 
zero zeta potential determined in this way agrees well with the values obtained in 
the literature. These experiments have resulted in a new and nondestructive method 
for determining the point of zero zeta potential of macroscopic semiconductor 
materials • 
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The effect of platinized powdered TiC >2 on the photoproduction of the halogens 
from aqueous solutions of their respective ions was investigated. It was found that 
platinized TiO^ powder was more effective in enhancing the desired reaction than 
unplatinized powder, and the reaction rates were higher at lower solution pH. Also, 
the production rate for the halogens is in the order I 2 > > CI 2 and is caused by 

the difference in the ability of the halide anions to compete with water for the 
photoproduced holes in the semiconductor. The demonstration of the photocataly tic 
production of chlorine and bromine from solutions of their ions in the presence of 
platinized powdered Ti0 2 suggests an inexpensive industrial method for producing 
these chemicals. 

Layered transition metal thiophosphate crystals were fabricated and tested for 
use as photoelectrodes in photoelectrochemical cells. Of the materials fabricated, 
only SnPS^, NiPS^, and FePS^ were stable, and these were stable only in acid solu- 
tions. Crystals of SnPS^ exhibited intrinsic n- and p-type photocurrent response. 

The iron-based and nickel-based sulfur compounds showed only p-type photocurrent- 
voltage characteristics. The materials had an indirect bandgap of about 2.1 eV and 
an onset potential near zero referenced to the saturated calomel electrode (SCE) at 
pH-2. These materials showed low conversion efficiencies when used as photoelec- 
trodes in liquid-junction photovoltaic cells with iodide and tri-iodide as a redox 
couple • 

The laser conversion efficiency of photoelectrochemical cells was explored using 
a technique of intercepting the light reflected from the window of the photoelectro- 
chemical cell with a photovoltaic cell to provide the external bias needed to obtain 
the optimum conditions for conversion. The technique was demonstrated using a-Fe 2 C >2 
as the photoelectrode and three silicon solar cells with laser radiation energy of 
2.4 eV. The silicon solar cells provided an external bias of 1 .0 V using only the 
reflected laser radiation. On the basis of this technique, an efficiency of 38 per- 
cent can be achieved for the conversion of laser radiation to storable chemical 
energy in the form of hydrogen and oxygen. This can be achieved by using SrTiO^ as a 
photoelectrode, zinc sulfide as the photovoltaic cell, and a laser photon energy of 
3.8 eV. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
October 29, 1982 
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SEMICONDUCTOR ENERGY-BAND MODEL 

The simple energy -band model of a semiconductor is outlined in this appendix and 
provides the background for our treatment of the semiconductor electrode. This out- 
line is included here to provide a concise statement of all the assumptions which 
underlie this work. Some of the details omitted in this appendix can be found in 
references 54 through 57. 

This model is based upon a simplified view of energy -band structure, which is 
applicable to all direct bandgap semiconducting materials. The actual band structure 
of a material used as a semiconducting electrode is approximated by two energy bands, 
each with a single extremum aligned as shown in figure 40. The energy -momentum rela- 
tion near an extremum at k = 0 can be approximated by the quadratic equation 


E(k) = E( 0) ± k 2 

2m 


where m* is the effective mass of particles occupying states in that band. This 
leads to an expression for the density of states in the conduction band 



( E > E c ) 
(E < E c ) 


and an expression for the density of states in the valence band 



(E > E v ) 
(E < E y ) 


The density of states is represented graphically in figure 40. 


(A1 ) 


(A2) 
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The equilibrium energy distribution of electrons about a Fermi energy E p is 
contained in the Fermi distribution function 


f (E) 



(A3) 


The concentration of electrons in the conduction band at thermodynamic equilibrium in 
the dark is given by the integral 


n 



p (E) f(E) dE 
r c 


(A4) 


and the concentration of holes in the valence band is given by the integral 


P 



p (E) [l - f(E) ] dE 


(A5) 


After a change in variable, equations (A4) and (A5) can be rewritten as 



(A6) 


(A7 ) 


F 1/2 (X) 



+ exp(x + s 


)] 


-1 


ds 
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A graph of this function appears in figure 41 • If the independent variable in the 
Fermi integral above is larger than 3, the integral can be approximated and evaluated 
in closed form as 


F 1/2 (X) 



expt-x 


s] ds 


jz 

2 


exp(-x) 


(A8) 


* Ep E„ E-p, E-- 

In this case the variable x is either ^ or — *-• If the Fermi energy 

is at least 3kT away from either the conduction band edge or the valence band edge, 
the concentration of electrons and holes can be approximated by the Boltzmann distri- 
bution as follows: 


n = N c 



with 


( A9 ) 


N 


C 



and 


P 




with 


(A10) 



At room temperature, the factor 3kT has a value smaller than 0.1 eV, so it is reason- 
able to use Boltzmann statistics in equations (A9) and (A10) to describe electron and 
hole concentrations at these temperatures. The product of the charge carrier con- 
centrations gives an expression that depends simply upon the energy gap, namely 


np 


W 



(All) 
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In an intrinsic semiconductor (chemically pure, single crystal), the only charge 
comes from thermally excited electrons and holes. The bulk of the semiconductor must 
be electrically neutral. Charge neutrality is specified by the condition 


n - p = 0 


which we express in terms of the intrinsic carrier concentration n^ 


by 


n = n ± = p 


(A12) 


Substituting the Boltzmann distributions for n and p, equations (A9) and (A10) 
into equation (A12), the Fermi energy for an intrinsic semiconducting material is 
related to the semiconductor parameters by 


E F,I 2^ + E c) 



(A13) 


or 


-r 2 ( E v + E C ) 


F , I 



(A14) 


This value of the intrinsic Fermi energy provides a reference level against which 
other energies can be conveniently measured. Its value is approximately half way 
between the energies of the valence and the conduction band edges so that it satis- 
fies the conditions required for the use of Boltzmann statistics. The intrinsic 
carrier concentration, defined by equation (A12), can be expressed in terms of the 
energy gap using equation (All) as follows: 


n i = exp 

In a semiconducting crystal that is not chemically pure or nonstoichiometric (a 
nonintrinsic semiconductor), the condition of charge neutrality in the bulk given by 
equation (A1 2) must be modified to reflect the statistical charge distribution asso- 
ciated with impurities or imperfections within the crystal. This modification 
depends upon structural details of the nonintrinsic semiconductor which are not often 
available. In this paragraph, a simple picture of donor- and acceptor -atom occupa- 
tion statistics that provides an approximate treatment of the problem of occupation 
statistics of impurity levels is described. This description is similar to that in 
references 57 through 60, but other more exact treatments can be found elsewhere 
(refs. 61 and 62). If the nonintrinsic effects are considered to be caused by a 
single localized one-electron state of energy communicating with the crystal by 



(A15) 
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electron transfer, the probability that the level is occupied by an electron or hole 
is given by 



where g is the number of degenerate states associated with the vacant level, and 
g* is the number of states associated with the occupied level. By representing the 
ratio of degeneracy in the form 



the probability of occupancy can be written in terms of a single parameter, the 
effective energy level E g , by the usual Fermi distribution function 


f 



1 



where 


(A16) 


* 

E 

s 


E + AE 
s 


E 

s 


+ kT In 



(A17) 


As shown in equation (A1 7) , the effect of the degeneracy of states at room tempera- 
ture is minimal. 

The semiconductor is assumed to have a concentration of acceptors with 

energy level E A just above the valence -band energy, and a concentration of donors 
N with energy level E^ just below the conduction-band energy. (See fig. 42.) 

The density of electrons occupying donor levels (un-ionized donors) n D and the 
density of holes occupying acceptor levels (un-ionized acceptors) p A obtained from 
equation (A16) are then expressed by the equations 


n^ = N 
D D 



(A18) 
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and 


-1 

(A19) 


The condition for charge neutrality in the bulk requires that the total number 4 

of negative charges, electrons plus ionized acceptors n + ” ^A be ec 2 ua ^ to the 

total number of positive charges, holes, and ionized donors p + - n D This condi- 
tion can be expressed by the equation r 



When the Fermi energy differs from either effective ionization energy (E^ or E D ) 
and is between them, complete ionization is expected, and in this case the condition 
for charge neutrality becomes 



n - p 


n d “ 


N 


(A21) 


Assuming that the electron and hole concentrations are governed by Boltzmann statis- 
tics (eqs. ( A9 ) and (A10)), the following relation involving the Fermi energy and the 
other characteristics of the semiconductor can be derived: 


sinh 



- n t 

D l 

2n. 

1 


( A22) 


where 



and 


6 
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By solving for the Fermi energy directly, we obtain for the case of an n-type 
semiconductor 


E 


F 


E 


F,I 


+ kT In 



and a p-type semiconductor 


E 


F 


E 


F,I 



is obtained. 


( A23 ) 


( A24) 
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TABLE I.- COMPARISON OF pH AT POINT OF ZERO ZETA POTENTIAL (pH ) DETERMINED 
FOR TWO-SOLUTION pH VALUES WITH VALUES CITED IN LITERATUrI FOR 
OXIDE-SEMICONDUCTOR ELECTRODES 


Electrode 

pH at 
25°C 

ref. RHE 3 
FB at 25°C 

Slope, b 
mV/deg 

n, c cm~^ 

N c , cm -3 

pH pzzp 

Average 

P H pzzp 

Literature 

P H pzzp 

Fe 2°3 

2 

0.5 

2.9 

5.1 X 10 17 

2.4 x 10 22 

9.8 

9.3 

8.6 


13 

.4 

2.7 



8.8 



SrTiO^ 

2 

-0.18 

2.25 

3.2 x 10 17 

6.8 x 10 19 

8.9 

9.3 

8.6 


13 

-.24 

2.4 



9.7 



Ti °2 

2 

0.11 

1 .15 

6.2 x 10 17 

3.0 x 10 21 

4.1 

5.0 

5.8 


13 

-.07 

1.5 



5.8 



wo 3 

2 

0.53 

0.1 

7.6 x 10 19 

9.8 x 10 18 

0.5 

0.5 

0.43 


Semiconductor parameter used in equation (8.2)* 

^Taken from figure 20. 

determined from Mott-Schottky plots. 


Ln 

cn 




TABLE II.- COMPARISON OF HALOGEN PRODUCTION FOR CATALYSTS AND pH VALUES 


Solution 

pH a 

Cat a lys t b 

Irradiation time, 
hours 

Amount produced, 
micromole 

1 .7 M KBr 

d 7 .0 

1 00 mg Pt-TiC >2 

1 

e 1 40 

1 .7 M NaCl 

d 7 .0 

100 mg Pt-Ti0 2 

1 

e 98 

2.0 M KI 

5.5 

None 

2 

1 .4 

2.0 M KI 

5.5 

150 mg AI 2 O 2 

2 

2.3 

2.0 M KI 

5.5 

1 50 mg Ti0 2 

2 

3.5 

2.0 M KI 

5.5 

150 mg Pt-Ti0 2 

2 

4.9 

2.0 M KI 

5.5 

150 mg Pt-Ti0 2 

4 

f 10.8 

2.0 M KI 

2.0 

150 mg Pt-Ti0 2 

2 

158 

3.0 M KBr 

5.5 

150 mg Ti 0 2 

40 

1 .3 

3.0 M KBr 

5.5 

150 mg Pt-Ti 02 

3 

.9 

3.0 M KBr 

5.5 

1 50 mg Pt-Ti0 2 

23 

6.2 

3.0 M KBr 

1.8 

150 mg Pt-TiC >2 

3 

1 1 

3.0 M NaCl 

5.5 

150 mg Pt-Ti0 2 

3 

0.9 

3.0 M NaCl 

5.5 

200 mg Pt-TiC^ 

24 

2.6 

3.0 M NaCl 

1 .2 

150 mg Ti0 2 

40 

4 

3.0 M NaCl 

1.2 

150 mg Pt-Ti0 2 

3 

13 

3.0 M NaCl 

1 .2 

150 mg Pt-Ti0 2 

20 

40 


a Adjusted by adding H 2 SO 4 or NaHCO^ * 

kpt-Ti0 2 prepared according to reference 34 with doped TiC > 2 anatase. 
Ti 02 powder was undoped anatase. 

C 1000-W high-pressure Hg lamp irradiated through glass window. 

^Adjusted by adding NaHCO^ • 

e Determined by initially introducing a known amount of As (III) to 
the cell. The As (III) also reacted, but at a lower rate than the halide 
as determined by blank experiments. 

^ Determined by continuously passing oxygen through the cell and 
through two external traps, each containing 35 ml of 2-M KI solution, 
and by titrating with As (III) both the I 2 in the irradiated cell and 
the I 2 trapped in the two external traps. 



TABLE III.- CRYSTAL GROWTH CONDITIONS, ONSET POTENTIALS 
BANDGAP, AND CONDUCTIVITY TYPE OF LAYERED 
TRANSITION-METAL THIOPHOSPHATES 



a Undertermined • 
b No photoresponse. 


TABLE IV.- MAXIMUM LASER-TO-CHEMICAL ENERGY CONVERSION EFFICIENCIES FOR 
SEMICONDUCTOR-BASED PHOTOELECTROLYTIC CELLS 


Semiconductor 


WOo 

Ti° 2 

SrTi0 3 

a-Fe 9 0 3 

BaTi0 3 

KTa0 3 

KTa 0.77 Nb 0.23°3 


V pB ref. RHE 

Bias, V 

Laser photon 
energy, eV 

Maximum 

efficiency, percent 

0.4 

4 

2.8 

6 

-.1 

1.5 to 2 

3.5 

33 

-.4 

1 .5 to 2 

3.8 

38 

.4 

2 

3.5 

15 

-.2 

1.5 to 2 

3.9 


-.3 

1 .5 

5.0 

15 

-.3 

1 .5 

5.0 

10 






’igure 1.- Schematic of a photoelectrochemical cell and current-voltaqe response 

for dark and illuminated conditions. 
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ELECTROLYTE 


1 SEMICONDUCTOR 


(a) Before charge redistribution. 


ELECTRON ENERGY 



(b) After charge redistribution. 

Figure 2.- Schematic of energy levels in an n-type semiconductor electrode 

and electrolyte. 
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Figure 3.- Energy-band diagram of an n-type semiconductor-electrolyte interface. 
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Figure 4.- Schematic of redox levels in electrolyte showing dissociation energy 

and surface states. 
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Figure 5.- Ionic charge distribution in the semiconductor-electrolyte interfacial 
region. Potential of this electrostatic field due to ionic and mobile charge 
distributions is shown in bottom half of sketch. An n-type semiconductor 
depletion region is assumed. 
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Figure 12.- Photocurrent as a function of incident laser power at 351 and 364 nm 
for a SrTiO^-based photoelectrochemical cell showing saturation. 
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Figure 13.- Ideal electrolytic voltage for water. 
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THERMAL PHOTOELECTROCHEMICAL CELL 


Figure 14.- Schematic of experimental setup of a thermal photoelectrochemical cell. 
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Figure 15.- Typical photocurrent-voltage curves for a SrTi0 3 ~based photoelectro lytic 
cell in 0.1 M NaOH for four temperatures. Potential is referenced to Hg/HgO 
electrode. 
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Figure 16,- Relative power conversion efficiency of an SrTiO^-based 
photoelectrolytic cell as a function of cell temperature. 
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Figure 17.- Photocurrent response at 1.5 V (referenced to Hq/HgO) to 
broadband radiation plotted aqainst cell temperature for SrTiO^ and 
Ti0 0 photoanodes • 
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Figure 18. 
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- Light-source spectrum and calculated thermal bandqap shift for SrTiO 

two temperatures. 



9.- Schematic of experimental setup to determine onset potential as a 

function of temperature. 
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Figure 20.- Temperature dependence of onset potential for Ti0 2 , SrTiO 

and Fe_0_. 




Figure 21.- Energy-level diaqram for a semiconductor-electrolyte junction in 
equilibrium with a reversible hydroqen electrode. 
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Figure 22.- Comparison of literature values of P H p ZZ p 
with values determined in this work as a function of 
electronegativity. Curve is taken from reference 24. 
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Figure 27.- Band edge of n-type SnPS 3 and p-type FePS 3 crystals 
in relation to redox levels of water, SCE, and the redox couple 
iodide and tri-iodide at pH-2. 
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Figure 28.- Spectral photocurrent response of n-type SnPS^ at 0.6 V referenced 

to SCE. 
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Figure 29.- Spectral photocurrent response of n-type SnPS^ at 0.6 
SCE. Linear dependence indicates an indirect bandqap, and the 
the bandgap energy. 
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Figure 30.- Spectral photocurrent response of p-type FePS^ at -1.5 V referenced to 
SCE. Straight-line dependence indicates an indirect bandqap, and the intercept 
gives the bandqap energy. 
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Figure 33.- Spectral response peak as a function of cell temperature for 

n-type SnPS Q . 
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Figure 34 •- Photocurrent density plotted against cell potential for n-type SnPS^ in a 
solution of 1 M H 2 S0 4' 0.025 M I^, and 1 M I - with a platinum counter electrode 
under white-light illumination. 
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Figure 40.- Conduction and valence band edges bent by deviations from charge 

neutrality near surface of semiconductor. Bendinq is depicted by plotting the 
energy of an electron against depth in semiconductor. Function <J> ( x ) is a 
normalized potential function used to describe shape of energy curve. 
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Figure 42.- Enerqy levels in a semiconductor showing conduction and valence bands, 
Fermi level, and location of donor and acceptor levels. 
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